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Abstract 

We show that, in the light-cone gauge, it is possible to derive in a 
very simple way the solution of the classical Yang-Mills equations for the 
collision between a nucleus and a proton. One important step of the 
calculation is the derivation of a formula that describes the propagation 
of a gluon in the background color held of the nucleus. This allows us to 
calculate observables in pA collisions in a more straightforward fashion 
than already proposed. We discuss also the comparison between light- 
cone gauge and covariant gauge in view of further investigations involving 
higher order corrections. 
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1 Introduction 

The study of semi-hard particle production in high energy hadronic interaction 
is dominated by interactions between partons having a small fraction x of the 
longitudinal momentum of the colliding nucleons. Since the phase-space density 
of such partons in the nucleon wave function is large, one expects that the 
physics of parton saturation nrm plays an important role in such studies. This 
saturation generally has the effect of reducing the number of produced particles 
compared to what one would have predicted on the basis of pQDC calculation 
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with parton densities that depend on x according to the linear BFKL mH] 
evolution equation. 

It was proposed by McLerran and Venugopalan that one could take 

advantage of this large phase-space density in order to describe the small x 
partons by a classical color field rather than as particles. More precisely, the 
McLerran-Venugopalan (MV) model proposes a dual description, in which the 
small X partons are described as a classical field and the large x partons act as 
color sources for the classical field. In their original model, they had in mind a 
large nucleus, for which there would be a large number of large x partons (at 
least iA where A is the atomic number of the nucleus, from just counting the 
valence quarks) and therefore they would produce a strong color source. This 
meant that one has to solve the full classical Yang-Mills equations in order to 
find the classical field. But that procedure, on the other hand, would properly 
incorporate the recombination interactions that are responsible for gluon satu¬ 
ration. In the MV model, the large x color sources are described by a statistical 
distribution, which they argued could be taken to be a gaussian for a large 
nucleus at moderately small x (see also |5] for a more modern perspective on 
that). 

Since then, this model has evolved into a full fledged effective theory, the 
so-called “Color Glass Condensate” (CGC) | 1()I11I12 |. It was soon recognized 
that the separation between what one calls large x and small x, inherent to the 
dual description of the MV model, is somewhat arbitrary, and that the gaussian 
nature of the distribution of sources would not survive upon changes of this sep¬ 
aration scale. This arbitrariness has been exploited to derive a renormalization 
group equation, the so-called JIMWLK equation |l:ill4IL^Il6llVll8ll9ll0lllll2| . 
that describes how the statistical distribution of color sources changes as one 
moves the boundary between large x and small x. This functional evolution 
equation can also be expressed as an infinite hierarchy of evolution equations 
for correlators m, and has a quite useful (and tremendously simpler) large Nc 
mean-field approximation EH, known as the Balitsky-Kovchegov equation. 

In high energy hadronic collisions, gluon production is dominated by the 
classical field approximation, and calculating it requires to solve the classi¬ 
cal Yang-Mills equations for two color sources moving at the speed of light 
in opposite directions. This is a problem that has been solved numerically in 
|22l28l24l25l26| for the boost-invariant case, and the stability of this solution 
against rapidity dependent perturbations has been investigated in |27l28j . In 
terms of analytical solutions, much less is known. The only situation for which 
gluon production has been calculated analytically is the case where one of the 
two sources is weak and can thus be treated at lowest order (this situation is 
often referred to as “proton-nucleus” collisions in the literature, but it can also 
be encountered at forward rapidities in the collision of two identical objects). 
This was done in a number of approaches The last two references 

provide the solution of Yang-Mills equation in this asymmetrical situation, in 
the Schwinger gauge {x'^A~ +x~A'^ = 0) and Lorenz gauge {df^A^^ = 0) respec¬ 
tively. More recently, Balitsky has proposed an expansion in commutators of 
Wilson lines, where at each order one treats the two projectiles symmetrically 
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Although the solution in the Lorenz gauge was fairly compact, it turned out 
that it displayed some oddities like the appearance of a Wilson line containing 
the coupling constant 5/2 instead of g (this of course disappeared at later stages 
from physical quantities). When applied to the production of quark-antiquark 
production in mni’ it also led to a contribution in which the vertex producing 
the qq pair was located inside the nucleus, which is quite counter-intuitive. 

In this paper, we derive the solution of the Yang-Mills equations in the light- 
cone gauge A+ = 0 (with the nucleus moving in the negative z direction), and 
we find a much simpler solution (not only the solution is simpler, but it is also 
much easier to obtain). In particular, the solution in the = 0 gauge presents 
none of the odd features encountered in the Lorenz gauge. Our central result, 
derived in section [21 is in fact a formula that tells how a color field propagates 
on top of the gauge field of the nucleus. The eqs. in fact contain all 
the information which is needed in order to derive the solution of Yang-Mills 
equations and compute gluon production, which we perform as a verification in 
section|21 Beyond the study of “proton-nucleus” collisions themselves, eqs. 
are an important building block for calculating higher order corrections in the 
weak source to the solution of Yang-Mills equations. And to a large extent, the 
complexity of this object finds its way into the solution of Yang-Mills equations. 
Therefore, it is important to determine this object, and to find a gauge in which 
it is particularly simple. For the sake of comparison, we derive in appendix 
m the analogue of eqs. m for the Lorenz gauge, and they are much more 
complicated, as one expects. 

2 Gluon propagation in a nucleus 

We assume that the nucleus is moving with a velocity very close to the speed 
of light, in the negative z direction. Thus, it can be described by the following 
density of color sources: 

S{x~^)pa{x±) . (1) 

Some intermediate calculations may require to regularize the delta function by 
giving it a small width. When this is necessary, we replace (5(a;“'') by (5e(a;“''), 
where Se(x~^) is a positive definite function normalized by dx~^ Se(x~^) = 1 
and whose support is [ 0 , e]. 

In this paper, we address the following question: knowing the gauge fields 
at x~^ = 0 , what are the gauge fields at = e, i.e. after having propagated 
through the nucleus? Finding the gluon propagator inside the nucleus amounts 
to solving this problem to linear order in the incoming gauge fields. In practice, 
we add to the gauge field Ag of the nucleus a small perturbation Al(, 

Af^ = A^ + A'^, ( 2 ) 

and we wish to find a linear relationship for this perturbation before and after 


3 


the region where the nucleus lives, such as (see figure : 

= e) = = 0) . 


( 3 ) 
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0 
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Figure 1: Gluon passing through a nucleus. The region shaded in gray is the 
region where the color source representing the nucleus is non-zero. 


We work in the light-cone gauge = 0. The covariant conservation of the 
color current reads^: 


d+J- + [D-,J+\ - [D\ /] = 0 . (4) 

This equation can be solved by 

J+ = J* = 0 , 

d+J- = 0 . (5) 

The first equation is allowed because all the color charges in our problem are 
moving in the negative z direction. The second equation means that the color 
current J~ associated to the nucleus is not affected by the incoming field^, and 
therefore this current J~ can be “hidden” once for all in the gauge field of 
the nucleus. 

The Yang-Mills equations in this gauge read 

d+{d^A>^) + ig[A\d+A] =0 

[D-,d+A-] - [D\F^-] = J- 

d+F-^ + [D -, d+A^] - [D^, FA] = 0 . (6) 

One can see that the first of these equations does not contain any time derivative 
{d~). Therefore, it can be seen as a constraint that relates the various field 
components at the same time. 

The gauge field of the nucleus alone is found by considering only the order 
zero in A^ in the above equations. One can readily check that the following is 
a solution: 

^0=0 , ^0 = = -9^{x+)^p(xi_) ■ (7) 

_ 

^In this and in the following equations, we do not write explicitly the color indices. 

^In other gauges, the incoming would induce a color precession of J~. 
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In order to find the the linear relationship between the perturbation 
before and after the region where the color sources of the nucleus live, we must 
linearize the Yang-Mills equations in One obtains 

- 5M1 = 0 

UA\ - 2ig[A-^,d+ A^] =0 

□Ar - 2ig[A^,d+A^] = 2i5[5Mo-, Al] = 2ig{d^A'^ ■ T)A\ . (8) 

The method for solving the second and third equations was explained in m 
The solution of the second equation inside the nucleus, i.e. for x~^ G [0, e], reads: 


A\{x'^,x ,x±) = U{x'^,0,x±)A\{0,x ,a:j_), 
where C/ is a Wilson line in the adjoint representation of the gauge group: 


U{x~^, 0, x±_) = 71|_ exp 


*5 


dz+ Ao-(z+,x^)r“ 


(9) 

( 10 ) 


At this point, we could find Aj” simply by solving the constraint equation. But 
as a verification, it is interesting to solve explicitly the third equation for Aj”, 
and at the end to verify that the fields A\ and Aj" do obey the constraint. The 
third equation leads to 

Ac(x+, a;“, a;j_) = U{x'^ ,Q,x±_)A^{Q,x~ ,x±_) 

+ig j dz+C/(x+,2:+,a;_L)(a*Ao (z+,a;_L) • r)C/(z+,0,a;_L)^Ai(0,x“,a;±) 

= [/(a;+,0, a:i)AC(0,a:“,a;i) -h {d'^U{x'^,Q,xx))^A\{Q,x~,xx) ■ (11) 

It is trivial to verify that the constraint that relates Aj" to A\ would have given 
the same answer. 

Therefore, if we denote U = C/(e, 0, x±) and if we don’t write explicitly the 
x~ and x± dependence, we have the following relations: 

A+(e) = 0 , 

Al(e) = [/A1(0) , 

Ar(e) = C/Ar(0) + (a*C/)^Al(0). (12) 

These equations are the light-cone gauge expression of the linear relation we 
were looking for. Analogous relations will be derived for the Lorenz gauge in 
the appendix fXl 


3 Gluon production in pA collisions 

3.1 Gauge field 

^From this linear relation, it is easy to calculate the gauge field that describes 
proton-nucleus collisions. In this case, the incoming field A^ is the field produced 
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by the current associated to the proton 


J+ = gS{x )pp{x±) . (13) 

For x~^ < 0, i.e. before the collision with the nucleus, the current J'*' remains 
constant, and we simply have, to linear order in the proton source Pp 

A+ = A^ = 0, 

□^1 = -e{x-)d^pp{xx). (14) 

Solving the latter equation gives the following field at = 0 : 

A+(0)=Ar(0)=0 , A\iO)=eix-)^ppixx). (15) 

The next step is to use the eqs. 03 in order to find the gauge field A^ imme¬ 
diately after the collision with the nucleus, i.e. at x'^ = e. 


A+(e)=0, 

A\{e) = U0{x~)-^pp{xx) , 

(e) = id"U)-^0{x-)-^pp{xj_) . (16) 

The final step is to find the gauge field for a;+ > e. The equation that governs 
the evolution of Al is 


□^1 - 2ig[A^,d+^,] = -|^J+ , (17) 

and the current J"*" at x'^ > e is modified by color precession in the nuclear field 

Ag : 

J~^{x~^ > e) = US{x~)pp{x±) . (18) 

This is a direct consequence of current conservation D~ = 0. We must 
solve eq. 03 with an initial condition at x'^ = e given by eq. 03. Using the 
techniques of m, we obtain 

^iix)= f dy-(fyj_Glix,y)2d+Uiy^)0{y-)-^pp{y^) 

Jy+=t Oy^ 

+ [ G°(a;,?/)6»(y-)a;(C7(yj^)pp(yj_)) , (19) 

Jy+>e 

where G°^{x,y) is the free retarded propagator obeying DxG°^{x,y) = S{x — y). 
The two terms in this solution are illustrated in figure |3 
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Figure 2: The two contributions to gluon production in pA collisions. Left: the 
proton source emits the gluon before the collision with the nucleus. Right: the 
proton source goes through the nucleus before emitting the gluon. The thick 
solid line represents the proton color current. 


3.2 Gluon production 

The amplitude for the production of a gluon of momentum p and polarization 
A is given by the Fourier transform of the amputated gauge field, contracted in 
the relevant polarization vector, 

Mx{p) = J d^x . (20) 

In light-cone gauge, the sum over the physical polarizations reads 

^ {p)ef^*{p) = -g^j , (21) 

A 

and for this reason we need only the transverse components of the gauge field 
when calculating gluon production. Note that since the Dalembertian 'OiA\(x) 
is bounded inside the nucleus, the region 0 < x'^ < e does not contribute in 
the limit e ^ 0. Therefore, we can take e ^ 0 and disregard the interior of the 
nucleus when we calculate the amplitude. We first obtain"^: 

DA\{x) = 2S{x~^)S{x~){U - l)^pp{x±) 

— 9{x~)9{—x^)d^pp{xj_) — 9{x~)9{x^)d^{Upp{x^)) , (22) 

and then the Fourier transform gives 
-p^A\{p) = -/A(,„t„„(p) 



P" 

kl 1 

1 (27r)2 

_{p+ + ie){p~ -1- ie) 

kf±\ 


Xpp(fcij_ ) [C/(fc2±) - (27r)25(fc2±)] . (23) 

^When we evaluate the Dalembertian of the gauge field in the region x"*" < 0, we must 
multiply eq. 113 by in order to restrict this term to the region x'^ < 0. Equivalently, 

we could simply subtract it from eq. GHJ in order not to overcount it in the region x'^ > 0. 
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In this equation, fc2_L = P± — ki± and ^proton (-P) is the Fourier transform of the 
gauge field of a proton alone, i.e. the Fourier transform of eq. m- It is easy 
to verify that this expression leads to the standard result for gluon production 
in proton-nucleus collisions^. 


4 Conclusions 

In this paper, we have obtained the solution of the classical Yang-Mills equations 
for proton-nucleus collisions in the light-cone gauge A~^ = 0. An important 
intermediate step is the “transfer matrix”, given in eqs. ca, which tells how 
a color field propagates through the nucleus on top of the color field of the 
nucleus. It turns out that this object takes an extremely simple form in the 
gauge A+ = 0, especially when compared to what one obtains in the (covariant) 
Lorenz gauge (see appendix IXIl . This transfer matrix is central in deriving the 
solution of Yang-Mills equations for pA collisions, and will be a crucial building 
block for calculating higher order corrections in the weak source to this solution. 
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A Gluon propagator in covariant gauge 

It is sometimes useful to have expressions for the propagation of a gluon field 
inside the nucleus in a covariant gauge = 0). We derive in this appendix 

the analogue of eqs. m for the Lorenz gauge. In this gauge, the color field of 
the nucleus is the same as its field in the A+ = 0 gauge: 

= A = 0 , ^0 = , (24) 

c'x 

thanks to the independence of the nuclear sources on x~. The Yang-Mills equa¬ 
tion that controls the evolution of Af is very simple and reads 

□A+-ig[Ao-,5+A+] , (25) 

and its solution can be written as 

_ At (e) = VAt (0) , (26) 

®In I29i:«)l . the same gauge was used but in a diagrammatic approach which deals with 
quark and gluon interactions with the nucleus, so that, the derivation of gluon production is 
more involved. 



where we use the same compact notations as in eq. 03 and where ^ is a Wilson 
line, 


V = T+ exp 


■ 9 
*2 , 


dz 


+ A~ 


^0a(^ 




(27) 


that differs from U only in the factor 1/2 in the exponential. 
The Yang-Mills equation for Al reads 


□A1 - 2tg[A^,d+A{] = ig{d^A^ ■ T)A+ - ig{A^ ■ T)9M+ . (28) 


The solution of this equation reads 

711(e) = UA\{0)+^^U<S>{^^A^ ■T)®V^At{0)-t^U<S>{A^ ■T)®-^{VAt{Q)) , 

(29) 

where we have used the following compact notation: 

Ui®A®U 2= f dx^ Ui{e,x~^,x±)A{x~^,x±)U2ix^,0,x±) , (30) 

Jo 

with Ui ^2 any pair of Wilson lines evaluated at the same transverse coordinate, 
and A the nuclear field or one of its derivatives. With this notation, one has 


d^U = igU®{d^A- ■T)®U , 

= i|Y® (5Mo -T)® Y , (31) 

and the slightly less obvious relation 

u -V = i^U ® (Ag • T) ® Y = i^V ^{Aq ■T)®U . (32) 

A proof of the latter formula was given in m Thanks to these relations, it is 
straightforward to simplify eq. (El into 


Al(e) = UA\{0) + d^ 




-U—A+iO) 


(33) 


One can see that all the convolutions between U’s and Y’s have disappeared 
from this expression. 

Let us consider finally the Yang-Mills equation for A/. This equation in¬ 
volves the current J~, and in the Lorenz gauge the field Af induces a color 
precession of J“, which can be seen as a correction of order 0{A^) to this 
component of the current. The equation that determines this correction is de¬ 
termined from the covariant current conservation, which can be satisfied by 


j+ = Ji = 0, 

5+Jf =^5[A+, Jg-] , (34) 


where Jg is the minus component of the nuclear current in the absence of any 
extra field. If we recall that J/" = —d^A^^ we can solve the last equation as 

Jr=*g(c?lAg-.T)Y^A+(0). (35) 
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The Yang-Mills equation for reads 

DA^ -^g[A^,^+A^] = J^ + {tgf{A^ -T^At 

-2ig{A^ ■ T)d-A+ - ig{d-A^ ■ T)At 
+2ig{d^A^ ■ r)Al + ig{A^ • T)aMl . (36) 

It is straightforward to solve this equation and write its solution as 


A^{e)=VA^{Q) 


+*|y 0 {dlA^ • T) ® V-^At{Q) + 0 (^0 • Tf 0 V-^At(Q) 

... , 1 -.1 


-igV®{A-^ ■T)®d- Y—A+(0) 


d+' 

1 


ilv®{d-A^ ■T)®V^A+{Q) 


+igV 0 (aMo . T) 0 [U^A\ (0) + {v^Af (0)) - U-^Af (0) 
-z|y 0 (Aq- • T) 0 a* [U^A\{0) + d^(v^At{0)) - U-^AtiO) 


(37) 

After some lengthy manipulations based on eqs. ED and lED, one can eliminate 
all the convolution products 0 and obtain 


Ai (e) — VA^ (0) -I- 


■ 1 a* 

d^U - V— 

a+ a+ 




a;(o) 
1 


\Av^ - (S-v)± + 8ir^] Af ( 0 ) ^ 


Eqs. (EHll . and (ISSl thus constitute the equivalent of eqs. CD in the Lorenz 
gauge. One can see in these expressions the advantage of working in the light- 
cone gauge: not only the formulas as much more compact, but in addition they 
do involve only the Wilson line U, and not V. As a side note, one can check 
that 

d^A^(e) = Vd^A^(O) , (39) 

i.e. the Lorenz gauge condition is satisfied at a;’*' = e provided that it was 
satisfied by the field incoming at a;+ = 0. 
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